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Abstract 

We suggest a disordered traffic flow model that captures many features of traffic flow. 
It is an extension of the Nagel-Schreckenberg (NaSch) stochastic cellular automata for 
single line vehicular traffic model. It incorporates random acceleration and deceleration 
terms that may be greater than one unit. Our model leads under its intrinsic dynam- 
ics, for high values of braking probability p, to a constant flow at intermediate densities 
without introducing any spatial inhomogeneities. For a system of fast drivers p — > 0, the 
model exhibits a density wave behavior that was observed in car following models with 
optimal velocity. The gap of the disordered model we present exhibits, for high values of 
p and random deceleration, at a critical density, a power law distribution which is a hall 
mark of a self organized criticality phenomena. 
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1 Introduction 



The investigation of traffic flow has attracted the interest of physicists already a long 
time ago. Different approaches have been proposed[l,2]. One can distinguish macroscopic 
and microscopic ones. In macroscopic models the traffic is viewed on the one hand as a 
compressible fluid formed by the vehicles and then analyzed using hydrodynamical fluid 
theories [3], on the other hand as a gas of interacting particles and then treated using 
kinetic theories of gases based on the Boltzmann equation[4]. 

In the car- following theories [5 and references therein], which are a typical example of 
microscopic approaches, individual vehicles are distinguished and the equation of motion, 
for each one, is the analogue of the Newton's equation. Many car-following models have 
been proposed depending on a sensitivity parameter and reaction time. But, they lead to 
unrealistic description of the behavior of free vehicles which have a large distance to the 
next vehicle ahead and suffer from serious problems in the low -density limit. To over- 
come these problems an optimal velocity V opt was introduced and a typical form [6] gives 
satisfactory results compared with empirical data. But, it has been found empirically 
that Vopt depends on the traffic state[7]. For a detailed study of all the above approaches 
and theories we refer the reader to the review articles Ref 5,8 and references therein. 

CA approach can be considered as a powerful tool, in statistical physics, to model local 
and nonlocal interactions [9]. Nowadays, the simulation of traffic using cellular automata 
approach[IO] stands out for its simplicity. It is a microscopic description of traffic flow 
governed by simple rules that each individual driver follows. 

The stochastic traffic CA model introduced by Nagel and Schreckenbergfll] (NaSch) is 
governed by simple rules. It simulates single-lane one way traffic and is able to reproduce 
the main features of traffic flow as backward moving shock waves and the so-called funda- 
mental diagram, J = J(p). In real traffic the system dynamics is very complex. However, 
the stochasticity introduced into the model takes into account of some events due to 
human driving: the maximum speed fluctuations, overreactions at braking and retarded 
acceleration. The NaSch model has been intensively studied using both analytical and 
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numerical methods[12]. Many extensions of this model have been established in order to 
understand the rich variety of physical phenomena exhibited by vehicular traffic [13, 14]. 
Some of these phenomena, observed in vehicular traffic under different circumstances, 
include transitions from one dynamical phase to another, criticality and self-organized 
criticality, metastability and hysteresis, phase-segregation, etc. We note that no one of 
these extensions is able to reproduce the main phenomena of real traffic flow all together. 

Since inhomogeneities have relevant effects on the systems dynamics, many kind of 
disorder were involved in traffic flow models. It was shown numerically that the NaSch 
model with a quenched random deceleration probabilities displays queueing of cars with 
a power law distribution of gaps between the cars at low densities[15]. The jammed phase 
behavior is similar to that observed in the standard NaSch model. The introduction of 
defects in NaSch model has a very high impact on the fundamental diagram and the dy- 
namics of the model[16,17]. One distinguishes between two kind of defects, the sitewise 
and the particlewise disorder. The later case produces phase separated stationary states, 
which consists of a large jam behind the slowest vehicle and a large gap in front of it at 
low densities. The single defect site induces a third phase, which is also a phase separat- 
ing, located between free flow regime and jamming phase where the flow is constant [17]. 
This intermediate phase was already observed in the asymmetric exclusion model, which 
corresponds to the NaSch model with v max = 1, for the sitewise disorder[18]. 

Since the CA approach may give, in a simple way, a good description of real traffic and 
captures the main phenomena observed by analyzing empirical data[7,19], we suggest an 
extension of the NaSch model that includes explicitly random incrementation of driver 
velocities. However, we think that it should be more realistic to take into account of non 
uniform acceleration and deceleration of vehicles as the drivers act differently. In the road 
you find careful drivers, which don't drive fast, as well as careless ones, which drive at 
their maximum speed. Depending on the distance ahead and the velocity of the vehicle 
the drivers may accelerate or decelerate more than one unit at each time step. Usually 
the competition between quenched randomness and dynamic fluctuations induces phase 
transitions between a disordered-dominated phase and fluctuation-dominated phase with 
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qualitatively distinct behaviors. Thus, we expect that a model with random acceleration 
should induce new dynamical states especially for intermediate densities since the fast 
drivers are stuck by the slowest ones. The model is defined in section 2 and depending 
either on random acceleration or deceleration we show numerically in section 3 that it 
presents some interesting phenomena which were observed separately in different varieties 
of traffic flow models[19-21,23]. In section 4 we study the gap distribution in order to 
show that the model exhibits a critical self organized behavior. The conclusion of our 
mean results is given in section 5. 

2 Definition of the model 

The disordered traffic flow model we present is a probabilistic CA where not only space 
and time are discret, but also the state variable of the vehicles. As in the NaSch model 
each cell can be empty or occupied by exactly one vehicle n and the state of each one is 
characterized by its velocity v n which can take one of the v max + 1 allowed integer values 
v = 0, 1, 2, .., v max . We denote the position and the velocity of the n-th vehicle by x n 
and v n respectively . Then, g n = x n+ \ — x n — I, where / is the vehicle length, is the gap 
between the n-th vehicle and vehicle n + 1 in front of it. At each time step t — > t + 1, the 
N vehicles arrange themselves on a finite lattice of length L following a parallel update 
according to the following rules: 
Step 1: Acceleration. 

If v n < v max , the velocity of the n-th vehicle is increased by a n sites, i.e. 



v n -> min(v n + a n ,v. 



max 



) 



(1) 



where a n = [p n g n ] + 1. 
Step 2: Deceleration. 

If (g n + 1) < v n , the velocity of the n-th vehicle is reduced to g n , i.e. 



v n min(v n ,g n ) 



(2) 
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Step 3: Randomization. 

If v n > 0, with probability p the velocity of the n-th vehicle is decreased randomly by d n 
sites, i.e. 

v n — > max(0, v n — d n ) with the probability p (3) 

where d n = [q n min([v max ], g n )\ + 1. 
Step 4-' Vehicle mouvement. 

Each vehicle moves forward according its new velocity v n obtained from the steps 1-3, i.e. 

x n > x n + v n (4) 

The symbol [A] denotes the integer part of A and the quenched p n and q n variables are 
randomly distributed in the interval [c, 1] according to the distribution laws: 

^) = (I^r(p- c ) B ( 5 ) 

and 

^)=( T ^t( 1 -^ (6) 
As the model we suggest is an extension of the NaSch model, the four steps quoted 
above are necessary to reproduce the basic features of real traffic [5, 13]. In step 1 the 
driver benefits from all the distance to the vehicle ahead and drives as fast as possible 
without crossing the maximum speed limit. He might move \p n gn] + 1 sites. In step 2, to 
avoid collision between vehicles, the driver reduces its velocity and adjusts it according 
to the distance to the vehicle ahead to be in security from any crash. The different 
behavioral patterns of the individual drivers as non deterministic acceleration as well 
as overreaction while slowing down are reflected in step 3. This step is crucial for the 
spontaneous formation of traffic jams. In real traffic the driver ability to drive more or 
less fast is a pertinent parameter in the dynamics of vehicles. The distributions given 
above reflect the fact that the 'careless' drivers drive as fast as possible. They correspond 
to higher values of p n and lowest values of q n . On the other hand the 'careful' drivers 
drive slowly since their corresponding probabilities p n and q n are the lowest respectively 
the highest. 
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3 Numerical investigations of the fundamental dia- 
gram 

Our numerical investigations are performed on a ring composed of L sites. The density 
of the vehicles is given by p = N/L where N is the number of vehicles. The flow < J > 
is expressed as p < v > where the mean velocity < v > of the vehicles is defined as 
the fraction of the sum of movable vehicle velocities among N vehicles. The length of the 
vehicles, /, will be taken as the unit of space coordinate. We use a parallel updating sheme 
since it takes into account the reaction-time and can lead to a chain of overreact ions. In 
what follow we will present the results obtained essentially from the distribution defined 
in eq(5,6) for n — 1. In order to analyze the effects of different random variables we will 
discuss 3 cases: 

3.1 deceleration effects 

We suppose that our vehicle drivers are all " careful" and they drive as slow as possible. 
They accelerate just by one unit at each time step while in the breaking step 3 they may 
reduce their velocity as maximum as possible. The fundamental diagram in this case, 
which corresponds to p n = for all vehicles, depends on the randomization parameter 
values p. However weak the value of p is, the flow, which usually presents a sharp max- 
imum, exhibits a smooth variation at its maximum value (Fig. la). This is due to the 
fact that the mean velocity of vehicles, v = v(p), is a decreasing function of p, even at 
low densities(Fig. lb). 

The existence of very slow vehicles in the system induces a jamming regions even at 
very low densities. The fast vehicles pile up behind the slowest ones. As a matter of fact, 
the space-time diagram presents, for intermediate values of p, a macroscopic high density 
region confined between relatively free flow ones(Fig. 2a). This results from the fact that 
the slowest vehicles, which are randomly distributed in the system, act like a blockage. 
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Figure 1: (a)The fundamental diagram of the model and (b) p-dependence of the velocity 



for v max = 5 in the case of random decelerations and uniform acceleration. 
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Thus, the flow gets a slower value than in the Nasch model. By increasing the global den- 
sity within the segregated phase the bulk densities of these regions remain constant, only 
their length changes. Consequently, the average flow is constant in the segregated phase, 
because the average density in the vicinity of the high density region does not depend on 
the global density. This intermediate phase, which is located for pi ow < p < phigh may 
be illustrated looking at the density profile which is plotted in fig 2b. For any particular 
choice of random q n variables according to the distribution law g(q)(eq. 6) one observes a 
separation into macroscopic high and low density regions. As the disorder average wasn't 
taken into account the density fluctuations observed in both regions result from the ran- 
dom decelerations of vehicles. They may also result from the interaction between vehicles 
that cause a formation of local jams (gaps) in the 'free' flow(congested)region. 

Thus, the average on the disordered q n variables leads to a quasi-constant flow in the 
segregated phase (Fig 3). As a result, one can distinguish, in the fundamental diagram, 
three different phases depending on the global density namely low density phase, high 
density phase and intermediate phase for p low < p < phigh- By increasing p in this inter- 
mediate phase, the flow presents a slow linear decrease. Its slope which is a decreasing 
function of p vanishes for some randomization parameter values. 

This phase separation located in our model was found in systems with ramps and systems 
with a stationary defect [8, 19]. For intermediate densities pi ow < p < pugh the flow is 
constant for these models. In this regime J(p) is limited either by the capacity of the 
ramp or the defect site. Even though the perturbations in these models are different, they 
all exhibit a plateau formation in the fundamental diagram as well as phase separation in 
the system independently of the nature of blockage. In the case of ramps it is the local 
increase of the density that decreases the flow locally. In the model with a stationary 
defect the increased slowdown parameter leads to a local decrease of the flow. In our 
model the "careless" drivers decelerate more than one unit and they act as a blockage in 
the system. Thus, the increased slowdown vehicles leads to a local decrease of the flow. 
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Figure 2: (a)The space-time diagram of the model and (b) the corresponding site- 
dependence of the density, in the case of random decelerations and uniform acceleration, 
for p = 0.3 and randomization parameter p = 0.8. 
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Figure 3: The fundamental diagram for p = 0.6 showing a plateau at intermediate densi- 
ties. 



Other variants of the model defined above may be obtained by modifying step2 or 
both step2 and step3 in the dynamics rules (eq. 2-3). If we assume that in the second 
step the vehicles may decelerate by reducing their velocities to be at a distance to the 
vehicle ahead more large than one unit such that v n — > min(v n , g n + 1 — d n ), the model 
presents a plateau formation in the fundamental diagram even for very low values of the 
randomization parameter p. The same result is obtained by changing both step2 and 
step3 of the model rules: 

step2: v n -> min(v n , g n + l- d n ) 

step3: v n — > max(0, v n — 1) with a probability p 

Since the deceleration effects in these model variants are emphasized, a plateau formation 
in the fundamental diagram occurs even at lower values of the randomization parameter 
p. 

3.2 Acceleleration effects 

Even though that the acceleration of vehicles is over estimated we think that it is 
worthwhile to study the effects of randomness of a such parameter. Those effects on the 
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fundamental diagram and on the dynamics of the model are studied by setting q n = and 
choosing randomly the variables p n for all vehicles. In this situation all the drivers are 
'careless' and drive as fast as possible. By performing numerical investigations we show 
that, for low values of p, the fundamental diagram(Fig. 4a) has quite similar form to that 
one obtained in the NaSch model. Thus, the distinction between these models should be 
formulated at microscopic level. 

As usual, the fundamental diagram depends on the randomization parameter p. For low 
values of p it presents two different regimes with respect to the density. At low density the 
system is in a free flow state and all the vehicles move, up to a certain density pi, at their 
maximum speed (Fig 4b). By increasing p this free flow region decreases considerably , i.e. 
Pi decreases, because the vehicles may frequently break and then reduce their speed more 
than one unit. This leads to the formation of spontaneous jams that are scattered in the 
space time diagram(Fig. 5a). By increasing the density more jams arise and are grouped 
together leading to a large strip of jammed region separated by free flow regions(Fig. 5b). 
Depending on the randomization parameter p the microscopic state of the system exhibits 
different structures. Indeed, for low values of p the system may selforganize into a density 
wave state(Fig. 6). For a such values of p most of the vehicles accelerate randomly and 
benefit from the distance of the vehicle ahead. In addition, the velocity of all vehicles 
in jammed region is limited by the slowest one. Then, if most of the vehicles located 
at the jammed region drive fast(slow)we have Jj n < J out (Ji n > J ou t), where Ji n (J ou t) 
denotes the number of entring(leaving) vehicles in the jammed region per unit of time. 
As the incoming vehicles pile up behind the jam we have p out < Pi n where p ut(pin) is the 
density of vehicles that leave(drive in) the jam. Consequently, the domain-wall between 
two stationary regions, free flow and congested region, propagates either in the opposite 
or the same direction of moving vehicles with the velocity v s h oc k, which depends on the 
values of Jj n and J out . From mass conservation we obtain 

_ jout Jin r Q \ 

V shock — — \o) 

Pout Pin 

As a result, the jammed regions wave from side to side. A such behavior was observed 
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Figure 4: (a) The fundamental diagram of the model and (b)p-dependence of the velocity 



for v max = 2 in the case of random accelerations and uniform deceleration. 
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in car following models with an optimal velocity that takes into account the character- 
istics of different vehicles[20]. In those models the different maximum velocities among 
vehicles are included in the optimal velocity function while in our CA model all vehicles 
keep the same v max and the characteristics of different drivers are induced in the acceler- 
ation terms. We note that the model presents regular oscillating behavior (Fig. 6a) for 
vanishing value of the randomization parameter, i.e. p — 0, while for finite values of p the 
breaking rule induces some irregularities. Detailed numerical investigations of the phase 
transition among the freely moving phase, the density wave phase and the homogeneous 
congested phase show that there is a critical value of the randomization parameter, p c 
above which the density wave does not appear [21]. If we denote by gf (gi) respectively Vf 
(vi) the average headway and velocity within (out of) the jam we may introduce an order 
parameter s = gf — gj or s = Vf — Vj and study the dynamics of the jammed region and 
the critical behavior of this inhomogeneous phase. 

3.3 Acceleration and deceleration effects 

In real traffic we found both "careful" and "careless" drivers. Then, both variables p n 
and q n should be chosen randomly. For very low values of randomization parameter p, 
most of vehicles drive fast. Consequently, at intermediate densities the steady state of the 
system reaches a density wave state with local 'defects' that are caused by spontaneous 
jams. This microscopic structure results from the interaction between the fast vehicles 
and the slowest ones. 

As in the case of slow drivers, i.e p n = 0, a region of constant flow occurs for high values 
of p(p > §). It becomes larger than in the previous case. Indeed, the faster drivers are 
blocked by the slowest ones and they pile up behind the jammed region. Consequently, the 
system exhibits, at intermediate densities, a high density band confined between relatively 
free flow regions. 
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Figure 5: The space-time diagram of the model, (a) p = 0.1 and (b)p = 0.2, in the case of 
random acceleration and uniform decelerations for v max = 5 and p — 0.8. 
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Figure 6: The space-time diagram of the model in the case of random decelerations and 
uniform acceleration for v max = 5, (a)p = 0.15 and p — 0, (b) p — 0.2 and p = 0.05 and 
(c)the spatio-temporal evolution of the gap, for p = 0.2 and p = 0.05. A density wave 
behavior is detected. 
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Figure 7: The gap distribution of the model, in the case of random decelerations and 
accelerations, for v max = 2, p = 0.05 and different values of p. 

4 The Gap distribution of the model 

Looking at the gap distribution over a wide range of densities in the congested flow 
region we see that it depends on the maximal velocity v max and the randomization param- 
eter p. In the case of high values of v max it exhibits, for low values of p, two maxima (Fig. 
7). The position of the maxima do not change considerably over a wide range of densities. 
The first maximum is assumed for the value zero. This means that the steady state of the 
system in the congested phase is of a densely packed queue and the two maxima of the 
gap distribution may be associated simply to a free flow phase with random distributed 
local inhomogeneities that are densely packed jams. We note that in a continuum limit of 
the NaSch model [22] the first maximum was observed at a non vanishing value and such 
a behavior was associated to a phase separation into congested and free flow regions. By 
decreasing the density, the first maximum vanishes and the gap distribution exhibits only 
one maximum that corresponds to the free flow regime. 

For high values of p located in a narrow interval, i.e. p ci < p < p C2 , the model displays, 
in the case of careful drivers, a power law gap distribution in the vicinity of some critical 
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density (Fig. 8). This behavior is a feature of the self organized criticality (SOC)[23]. The 
distribution doesn't change considerably with the system size and finite size data collapse 
of the form 

P(g) = g- T f(g/Ln (9) 

Using a finite size analysis it's easy to see that P(g) = L^ 13 f(g/ L u ) with /3/v — r. The 
critical exponents r and v vary continuously with the randomization parameter p and the 
maximum velocity v max . The SOC behavior results from the fact that the fast vehicles 
built behind the slowest ones forming platoons separated by large gaps. The power law 
decay of the gap distribution was illustrated analytically for the asymmetric exclusion 
model with random rates using some statistical estimates [24]. The numerical estimation 
of the critical densities and the critical exponent r of the gap distribution of our model are 
represented in table 1 for v max = 2 and v max = 5. We note that the SOC behavior was also 
observed in a stochastic traffic model with random deceleration probabilities [15]. It was 
pointed out that the NaSch model with random randomization parameter probabilities 
self organizes into a stable queueing phase at low densities and has a power law gap distri- 
bution. This model belongs to the same universality class of the BFL model [25] while our 
model belongs to another universality class as the critical exponents r are rather different. 



^max 


p 


Pc 


T 


2 


0.7 


0.19 


0.755 




0.8 


0.18 


0.812 


5 


0.75 


0.1 


0.915 




0.8 


0.11 


0.972 



Table 1: The estimations of the critical exponents r of gap distributions at critical den- 
sities p c for different v max and different values of the randomization parameter p. 
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Figure 8: The log-log plot of the gap distribution for v max — 2, p c — 0.1 and p = 0.75, in 
the case of random accelerations and uniform deceleration, and for different system sizes 
L ranging from L = 860 to L — 3440. The finite size scaling analysis gives a good fit of 
our data with the exponents r = 0.903 and v = 0.083. The inset shows the finite size fit 
with = 0.075. 

In the case of 'careful' drivers, all vehicles accelerate just by one unit at each time 
step but they may decelerate by reducing their velocity as maximum as possible. The 
probability that one driver decelerates as maximum as possible is realized for 1 — 8 < q n < 
1, where 5 « 1 and is calculated from the distribution given in eq (6): p m ~ 5 2 /2. At low 
densities the system evolves according to its own dynamics and is however interrupted by 
relatively small perturbations at vanishingly small rates. Consequently, the system gets 
enough time to organize itself and relax back to its corresponding steady state before it is 
perturbed again. This effectively separates completely the time scales for perturbing the 
system and its response. This infinite time scale separation is required by SOC behavior. 
The system is then self similar and the scale invariance is due to the fact that the jams 
are fractal in the sens that there are smaller sub-jams inside large jams. 
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5 Conclusion 



The dynamics of the disordered traffic flow, which is based on the NaSch model, pre- 
sented in this paper bears considerable structures and many of interesting features of 
traffic flow. Random acceleration and deceleration terms were introduced in the CA 
rules. Within numerical investigations we have shown that the model presents under its 
own dynamics many interesting phenomena that were observed separately in different va- 
rieties of traffic flow models namely the segregated phase, the density wave state and the 
SOC behavior. Depending on the randomization parameter p a rich variety of fundamen- 
tal diagrams is obtained. For intermediate values of p and low densities the fundamental 
diagram may trigger transitions from free to congested flow under its own dynamics with- 
out introducing any external perturbation like defects or on and off-ramps[7,13]. We may 
assume that at intermediate densities the slow vehicles act like a defect site leading to a 
blockage. As a matter of fact an intermediate phase, where the current is quasi-constant, 
is inserted between free flow and jamming phases. The analysis of the microscopic struc- 
ture of the space time diagram for low values of p shows that the system may organize 
itself into a wave density phase where the jams (gaps) oscillate. Such behavior was pointed 
out by Nagatani in the optimal velocity model [20]. The gap distribution of our model 
in the case of slow drivers exhibits, for some values of randomization parameter p and 
at low densities, a power law behavior in the vinicity of a critical density signaling the 
existence of SOC phenomenon. The critical exponents are rather different from those 
obtained in a disordered traffic flow model where the quenched disorder was incorporated 
in the random deceleration probabilities [18]. The model we suggest doesn't belong to the 
same universality class of the BFL model as the former model does. Finally, we note that 
the disorder introduced in our model is the essential ingredient to observe all the traffic 
features mentioned above since they don't occur in the case of uniform acceleration and 
deceleration[26],i.e p n = est and q n = est. 

The model we present in this paper presents many interesting features of traffic flow 
and may be developed in order to capture other pertinent phenomena observed in traffic 
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flow as the hysteresis and the synchronized state. 
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